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摘要 : 在 天 文 观测 中 ， 长 缝 光谱 仪 是 一 种 常用 的 光谱 观测 设备 。 国 家 天 文 台 兴隆 观测 基 
地 2.16m 望远镜 搭载 的 北京 暗 天 体 摄 谱 仪 (Beijing Faint Object Spectrograph and Camera, 
BFOSC) 的 长 颖 光谱 模式 应 用 广泛 ， 但 没有 配备 专用 的 数据 处 理 软件 ， 通 常用 户 观测 后 采用 
IRAF( Image Reduction and Analysis Facility) 软件 利用 交互 方式 处 理 。 介 绍 了 一 款 针 对 BFOSC 
长 颖 光谱 的 自动 化 数据 处 理 软件 ， 数 据 处 理 过 程 包 括 本 底 ( bias) 和平 场 (flat) 改正 、 背 景 改 


正 、 抽 取 光 谱 和 波长 定 标 等 ， 展 示 了 对 实测 数据 进行 处 理 的 结果 ， 并 与 RAF 软件 处 理 的 结 


果 进 行 了 对 比 。 该 软件 在 完成 所 有 光谱 处 理 必需 步骤 的 基础 上 实现 了 自动 化 


化 抽 谱 方法 ， 能 够 自动 扣除 大 部 分 宇宙 线 ， 保 证 了 高 精度 的 输出 结果 。 
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二 /一 


运行 ， 并 采用 优 


我 国 自主 研制 的 第 1 台 2 m 级 望远镜 一 国家 天 文 台 兴 隆基 地 的 2. 16 m 望远镜 ， 配 备 了 北京 暗 
天 体 摄 谱 仪 "1 。 这 人 台 仪 器 具备 5 种 不 同 的 观测 模式 ， 其 中 长 颖 光谱 是 应 用 最 为 广泛 的 一 种 模式 ， 可 以 
在 不 同 的 颖 宽 、 光 栅 和 波光 片 组 合 之 间 进 行 切换 。2016 年 更 换 CCD 相机 后 最 新 的 色散 、 波 长 范围 等 
参数 如 表 1553] 。 据 公开 资料 统计 ，2021 年 兴隆 观测 基地 为 2. 16 m 望远镜 用 户 分 配 的 252 个 观测 夜 


中 ， 使 用 BFOSC 的 有 165 夜 ， 占 比 超过 65% 。 
BFOSC 测 光 模式 下 的 数据 可 以 使 用 文 [4] 开 
发 的 程序 处 理 。 但 长 颖 光谱 部 分 尚 没有 专门 的 数 


#1 BFOSC 不 同 光栅 及 基本 性 能 参数 
Table 1 Grating parameters of BFOSC 


据 处 理 软件 ， 观 测 者 一 般 自 行使 用 IRAFI5 5 等 Name Rec. Lin. Disp. Disp. Wav. Range 
软件 利用 交互 方式 对 原始 光谱 图 像 进行 处 理 ， 获 /nm mm) /(nm. pix’) /nm 
得 一 维 谱 。 但 这 类 软件 并 不 是 针对 BFOSC 特点 G3 13.9 0. 208 330-660 
WIW, RAET KAE IERE AY ETT OEE, G4 19. 8 0. 297 360-870 
且 美 国 国立 光学 天 文 台 (National Optical Astronomy C5 19.9 0. 298 520-1 000 
Observatory, NOAO) 已 经 宣布 自 2013 年 起 不 再 C6 8.8 0. 132 330-545 
对 IRAF 进行 官方 支持 。 G7 9.5 0. 142 387-676 
ae ee a a a G8 8.0 0. 120 580-828 
详 处 理 软件 ， 实 现 了 自动 化 处 理光 详 数 据 。 该 程 a a6 — hag 
序 可 以 读 取 用 户 指定 的 光谱 仪 原始 数据 ， 自 动 进 
行 图 像 合 并 、 波 长 定 标 、 目 标 源 光 谱 抽取 等 步 。 “| 
又 ， 在 关键 步骤 中 保存 中 间 结果 并 进行 可 视 化 。 ”2 ee 
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程序 引入 了 优化 抽 谱 方法 ， 上 自动 屏 珊 了 宇宙 线 ， 提 高 了 最 终 数 据 产 品 的 精度 。 用 户 可 以 使 用 该 程序 快 
速 实 现 批量 数据 处 理 ， 大 大 减轻 了 用 户 的 工作 量 ， 有 效 提升 了 工作 效率 。 


1 系统 整体 框架 


1.1 环境 配置 

本 程序 基于 Python3 编程 语言 编写 ， 依 赖 
一 些 开 源 和 免费 的 第 3 方 软件 包 ， 包括 Astropy 
(用 于 FITS 文件 的 读 取 与 存储 ) ，Numpy，Scipy 
(用 于 科学 计算 ) 和 Matplotlib( 可视化 ) 。 本 软件 
在 www. zenodo.org/record/7866030#. ZEi9uM5ByUk 
以 开源 形式 发 布 。 
1.2 程序 结构 与 流程 

程序 主要 包括 数据 读 取 、 数 据 合 并 和 光谱 
抽取 3 部 分 。 其 中 数据 读 取 包括 扫描 生成 观测 
日 志 、 读 取 各 个 步骤 数据 。 数 据 合 并 包括 本 底 
合并 、 平 场合 并 与 改正 以 及 定 标 灯 红 蓝 端 合并 。 
光谱 抽取 包括 本 底 改 正 、 平 场 改 正 、 目 标 星 定 
位 、 岁 像 畸 变 改正 、 背 景 改 正 、 波 长 定 标 和 一 


Raw Observation 
images log 


Scientific 
objects 


Image 
combination 


Master 
bias 
Image 
combination 
Master Flat field correction | 
flat field I I 


| Curvature correction | 


Bias correction | 


维 谱 抽取 等 步 又 。 具 体 流程 如 图 1。 用 户 只 需 
要 输入 原 始 数据 文件 所 在 的 路 径 j 程序 A 动 执 Background subtraction 
行 必要 的 步骤 并 存储 各 步骤 产生 的 中 间 结 果 ， aan | ce 
以 便 后 续 使 用 对 关键 步骤 进行 可 视 化 并 得 到 extraction location 
最 终 数据 产品 。 I I 


Optimal 
extraction 


其 中 本 底 (bias) 是 指 没有 光子 输入 时 ，CCD 
像素 在 通电 时 ( 零 秒 曝 ) 的 读 出 值 。 平 场 (flat) 是 


Wavelength 
calibration 


8 输入 光源 为 均匀 光 的 情况 下 CCD 获得 的 图 
像 ， 反 映 了 CCD 的 各 个 像素 对 光子 的 响应 差 scientific objects 
异 。 定 标 灯 是 指使 用 空心 阴极 射线 灯 的 光谱 ， 

包含 波长 已 知 的 发 射 谱 线 ， 用 于 标定 CCD RR ae eee 

上 波长 与 像素 之 间 的 对 应 关系 。 科 学 目标 是 光谱 Be eee 


仪 的 使 用 者 拍摄 的 天 体 目标 ， 例 如 彗星 、 恒 星 和 星系 等 。 
1.3 程序 使 用 

本 文 介绍 的 程序 是 针对 BFOSC 的 长 颖 光谱 设计 ， 用 于 将 光谱 仪 获取 的 原始 观测 数据 (通常 为 
FITS 格式 的 光谱 图 像 ) 转 换 为 可 供 科 学 分 析 的 一 维 光 谱 ， 即 目标 天 体 的 波长 与 流量 的 对 应 关系 ， 最 终 
输出 数据 以 FITS 格式 保存 。 用 户 需 要 在 操作 系统 的 终端 运行 本 程序 ， 并 输入 原始 数据 文件 和 观测 日 
志 所 在 的 路 径 ， 程 序 自动 执行 必要 的 步 又 ， 存 储 中 间 结 果 并 进行 可 视 化 ， 得 到 最 终 数据 产品 ， 即 目标 
天 体 的 一 维 谱 。 这 些 结果 和 中 间 过 程 文件 都 存储 在 指定 的 计算 机 目录 中 。 此 外 ,用户 也 可 以 更 改 运行 
模式 ， 手 劲 对 目标 星 进行 标记 ， 程 序 会 根据 标记 位 置 寻找 目标 星 准确 轮廓 。 


2 实际 处 理 问 题 与 解决 方法 


2.1 扫描 生成 观测 日 志 
2. 16m 望远镜 的 观测 者 在 观测 同时 经 常会 记录 电子 版 的 观测 日 志 , 但 是 为 了 避免 格式 不 一 致 ， 
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程序 会 扫描 原始 数据 所 在 的 计算 机 目录 ， 读 取 其 中 FITS 格式 的 文件 ， 根 据 每 个 文件 的 头 部 (Header) 
言 息 ， 结 合用 户 记录 的 观测 日 志 ， 重 新 生成 可 以 机 读 的 格式 化 观测 日 志 ， 包 含 文件 名 、 图 像 类 型 、 观 
测 目标 源 的 名 称 、 曝 光 时 间 、 观 测 开始 时 间 等 信息 。 
2.2 本 底 改 正 

通常 情况 下 ， 观 测 者 在 观测 开始 前 或 结束 后 在 CCD 相机 快门 关闭 的 情况 下 进行 多 次 零 秒 曝光 ， 
获得 一 组 本 底 图 像 ， 通 过 图 像 合并 获得 本 底 的 数值 并 抑制 CCD 读 出 噪声 带 来 的 误差 。 这 个 步骤 提供 
了 两 种 不 同 的 合并 模式 ， 即 均值 模式 和 中 值 模式 ， 通 常情 况 下 默认 为 均值 模式 。 本 底 合 并 后 生成 本 底 
图 像 ， 并 将 合并 后 的 本 底 图 像 以 FITS 格式 单独 存储 为 一 个 文件 。 之 后 对 需要 操作 的 数据 直接 扣除 保 
存 的 本 底 合 并 数据 即 完成 本 底 改正 。 
2.3 目标 源 定位 

本 文 介绍 的 程序 采用 以 下 方法 确定 目标 天 体 的 光谱 在 CCD 图 像 上 的 位 置 : 从 光谱 图 像 的 左 侧 开 
台 ， 沿 x 轴 方 向 每 隔 50 个 像素 取 图 像 的 截面 最 高 处 作为 初始 位 置 ， 取 上 下 各 20 像素 范围 内 的 数据 ， 
用 高 斯 函数 进行 拟 合 ， 获 取 纵 向 (y 方向) 的 中 心 位 置 与 半 高 全 宽 ， 并 将 中 心 位 置 当 作 新 的 初始 位 置 。 
重复 上 述 过 程 ， 直 到 拟 合 参 数 收敛 。 接 下 来 对 所 有 y 方向 的 中 心 位 置 与 横 坐 标 (x 值 ) 的 关系 进行 三 阶 
多 项 式 函 数 拟 合 ， 作 为 目标 源 光谱 的 中 心 位 置 曲线 。 光 谱 图 像 如 图 2。 此 外 程序 还 为 用 户 提 供 了 手动 
选取 目标 源 的 模式 ， 在 该 模式 下 用 户 可 以 在 弹出 的 图 形 界面 上 用 鼠标 标注 目标 源 的 大 致 位 置 ， 程 序 在 
其 标记 的 位 置 附近 按 上 述 流 程 寻找 目标 天 体 的 光谱 。 


Trace for 202304180046 SPECLTARGET SN2023eoc slit18 385SLP_G4 (SN2023eoc) 
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图 2 (a) 中 心 位 置 与 半 高 全 宽 图 ; (b) 平 均 轮 廊 图 ; (c) 真 实 目标 源 位 置 图 ; (d) 轮 廊 中 心 拟 合 图 
Fig.2 (a) The central positions and full-width half-maximum (FWHM) of the spectral profiles, (b) the average profile is plotted 


in the lower left panel with black curve. (c) the target image, and (d) the positions and FWHMs of the center fittings 
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2.4 平声 改正 

由 于 CCD 的 各 个 像素 对 光子 的 响应 略 有 差异 ， 即 便 输入 光源 为 均匀 光 ，CCD 获得 的 图 像 也 是 不 
均匀 的 。 在 拍摄 过 程 中 一 般 进 行 多 次 曝光 ， 将 获得 的 多 幅 平 场 图 像 进行 琶 加 ， 以 获得 较 高 的 信 噪 比 并 
抑制 读 出 噪声 。 然 后 使 用 三 次 样 条 插值 补 齐 数据 中 的 坏 像 元 数据 ， 之 后 使 用 Savitzky-Golay 滤波 过 
Fe? 获得 主 平 场 图 像 。 每 一 幅 目标 天 体 的 光谱 图 像 需 要 除 以 主 平 场 图 像 进 行 平 场 改 正 。 
2.5 视 场 光学 畸变 改正 

由 于 CCD 相机 的 光学 畸变 ，CCD 图 像 上 狭 缝 所 成 的 像 存 在 一 定 程 度 的 弯曲 。 在 波长 定 标 灯 的 光 
谱 图 像 上 选取 沿 着 y 轴 中 间 位 置 的 一 维 谱 作为 基准 光谱 ， 在 y 轴 的 上 下 两 个 方向 ， 每 间隔 200 个 像素 
进行 一 次 同样 的 一 维 谱 提取 操作 ， 将 得 到 的 光谱 与 基准 光谱 计算 交叉 相关 函数 


LEO -用 (es =) 
ws ， 1 
/SD_SD, as 
其 中 ,x ly 分 别 为 待 改 正 的 光谱 与 基准 光谱 ，n 为 光谱 取 点 个 数 ; 上 为 相对 像素 差 值 ; SD. 和 SD, 分 
别 为 两 光谱 的 标准 差 。 交 又 相 关 函 数 最 大 值 对 应 的 像素 之 差 即 为 实际 的 像素 偏 移 量 。 根 据 各 行 的 像素 
差 获得 改正 函数 ， 并 生成 改正 图 像 ， 如 图 3， 不 同 颜色 的 线 代表 方向 不 同位 置 的 光谱 。 经 过 视 场 畸 
变 改 正 后 能 有 效 提高 一 维 光 谱 的 波长 定 标 精度 ， 吕 免 系统 误差 。 
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图 3 (a) RAE GH FeAr 定 标 灯 谱 的 一 维 谱 图 像 ; (b) R AER ao — eh AR 
Fig.3 (a) The FeAr spectra at different positions along the slit direction after the 


distortion correction; (b) before the distortion correction 


2.6 背景 改正 
实际 观测 过 程 中 ， 我 们 获得 的 光谱 不 仅 来 自 于 目标 源 ( 图 4(a) )， 也 有 天 光 背 景 以 及 仪器 内 部 的 
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杂 散 光 。 目 标 星 的 光谱 也 包含 这 些 成 分 ， 因 此 为 了 获取 目标 源 的 真实 流量 ,我 们 需要 对 图 像 数 据 进行 
背景 改正 。 

程序 一 般 会 控制 目标 星 使 其 落 在 沿 y 轴 中 间 附 近 的 位 置 。 程 序 自 动 选 取 CCD 图 像 上 距离 目标 星 
EE, CCD 像素 一 致 性 比较 好 的 行 (例如 图 4(b) 中 第 340~500 行 以 及 1600~1820 行 ) 作 为 测定 背 
景 光 大 小 的 区 域 。 在 选取 过 程 中 ,程序 会 计算 该 区 域 每 一 行 沿 着 x 轴 ( 即 色散 方向 ) 的 平均 值 和 标准 
差 ， 对 于 强度 超过 均值 以 上 3 倍 标准 差 的 行 ， 可 能 包含 目标 源 附近 其 他 天 体 的 光谱 ， 这 种 情况 在 密集 
星 场 中 尤为 常见 。 程 序 在 计算 背景 值 时 对 这 些 行进 行 了 屏蔽 ， 以 免 获得 比 真 实 值 显著 偏 高 的 背景 流 
量 。 目 标 源 光谱 附近 的 像素 值 减 去 背景 光谱 ， 即 为 目标 天 体 的 光谱 。 


202304180046_SPECLTARGET_SN2023eoc_slit18_385LP_G4 (SN2023eoc) 
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图 4 (a) 目 标 源 天 体 的 原始 光谱 图 像 ; (b) 本 例 中 测定 天 空 背 景 采 用 的 像素 区 域 
Fig.4 (a) The original 2D image of the target source; (b) the CCD image areas that are 


used to determine the background level 


2.7 波长 定 标 

波长 定 标的 目的 是 在 CCD 图 像 上 建立 像素 位 置 与 波长 的 映射 关系 。 本 文 介绍 的 程序 包含 自动 对 波 
长 进行 标定 ， 具 体 处 理 流程 如 下 : (1) 读 取 该 模式 下 提供 的 已 经 定 标 完成 的 模板 光谱 数据 、 模 板 特征 
发 射线 波长 与 像素 位 置 对 应 数据 ， 其 次 对 要 处 理 的 光谱 数据 与 模板 光谱 数据 进行 交叉 相关 函数 计算 ， 
获得 两 个 光谱 数据 位 置 的 偏 移 量 ， 将 读 和 人 的 模板 像素 位 置 数 据 增加 偏 移 量 获得 新 的 像素 位 置 数 据 。 
(2) 依 次 取 新 的 像素 位 置 左右 各 10 个 像素 的 一 维 光谱 ， 用 广义 高 斯 函数 4exp[ -(a-c)/a] +B 进行 拟 
合 ， 获 得 峰值 位 置 对 应 的 像素 数值 。 结 果 如 图 5 和 图 6， 其 中 每 幅 图 中 的 圆 点 代表 一 维 光 谱 数 据 ， 实 
线 代 表 拟 合 曲 线 ， 虚 线 代 表 拟 合 得 到 的 谱 线 中 心 位 置 。(3) 用 五 阶 多 项 式 函 数 进行 拟 合 ， 获 得 波长 与 
像素 位 置 的 对 应 关系 函数 ， 并 保存 全 部 波长 与 像素 位 置 数据 组 ， 夯 出 波长 定 标 函数 图 像 ， 如 图 7。 

此 外 ， 由 于 BFOSC 采用 的 定 标 灯 本 身 的 特性 ， 其 红 端 存在 大 量 明亮 的 Ar 元 素 的 发 射线 ， 而 蓝 端 
以 Fe 线 为 主 ， 它 们 的 强度 差异 较 大 ， 如 果 曝 光 时 间 过 长 会 导致 红 端 饱和 ， 而 上 曝光 时 间 过 短 会 导致 蓝 
端 谱 线 没有 足够 的 强度 。 因 此 ， 我 们 通常 拍摄 一 次 长 曝光 和 一 次 短 曝光 ， 分 别 作为 蓝 端 和 红 端 的 波长 
定 标 谱 。 在 本 文 的 例子 中 ， 红 蓝 两 端的 波长 定 标 图 像 曝光 时 间 分 别 为 30 s 和 300 s。 因 此 ， 在 进行 定 
标 之 前 本 方法 对 两 组 光谱 数据 进行 截取 之 后 合并 处 理 生成 一 个 新 的 定 标 光谱 图 像 ， 处 理 过 程 中 计算 交 
叉 相关 函数 时 对 红 端 和 蓝 端 分 别 求解 计算 。 
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AS 用 于 波长 证 认 的 FeAr 谱 线 峰 值 位 置 拟 合 
Fig.5 Peak positions of the identified FeAr lines 
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图 7 (a) 多 项 式 拟 合 得 到 的 波长 (y) 与 像素 (*) 的 关系 ， 蓝 色 点 代表 定 标 谱 线 ; (DURE 


Fig.7 (a) The relation of fitted wavelength v. s. pixel numbers, where blue dots represent the fitting emission lines; (b) the residuals 


优化 抽 谱 方法 采用 加 权 求 和 ， 保 证 具有 足够 流量 的 情况 下 包含 尽量 少 的 噪声 像素 ， 对 暗 弱 源 尤其 
有 效 ， 能 显著 提高 信 噪 比 ” 。 优 化 抽 谱 过 程 主要 分 为 两 步 : (1) 获 得 目标 星 的 精确 空间 轮廓 ， 确 定 每 
个 像素 的 加 权 值 ; (2) 根 据 空间 轮廓 的 形状 ， 排 除 明显 高 于 轮廓 的 像素 点 ， 即 扣除 宇宙 线 。 在 我 们 开 
发 的 光谱 数据 人 处理 程序 中 ,第 1 步 首先 默认 空间 轮廓 宽度 为 36 个 像素 ,轮廓 方向 每 个 像素 具有 相同 
的 权重 值 ， 然 后 计算 光谱 数据 方差 ， 根 据 方差 改正 加 权 值 与 空间 轮廓 ， 重 复 此 过 程 和 迭代 2~3 次 ， 最 
终 获 得 精确 的 空间 轮廓 。 第 2 步 ， 先 根据 上 一 步 得 到 的 轮廓 提取 科学 目标 的 一 维 谱 ， 并 与 普通 抽 谱 方 
法 得 到 的 一 维 谱 计算 差 值 ， 对 其 中 大 于 5 倍 标准 侦 差 的 像素 点 进行 屏蔽 ， 我 们 认为 这 些 像素 是 受到 宇 
宙 线 照射 的 像素 。 重 复 上 述 步骤 直至 收敛 ， 至 此 我 们 获得 了 屏蔽 宇宙 线 后 的 科学 目标 一 维 谱 。 
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3 结果 比较 


为 了 验证 程序 的 可 靠 性 ， 我 们 使 用 IRAF 对 超新星 SN2023eoc 的 同一 组 光谱 数据 进行 了 处 理 并 绘 
制 了 一 维 谱 图 像 ， 如 图 8， 其 中 ， 黄 色 线 是 本 软件 自动 化 处 理 获得 的 结果 ， 蓝 色 线 是 IRAF 处 理 结 
由 图 8 可 以 看 出 ， 两 种 处 理 方 法 结果 基本 一 致 ， 并 且 本 软件 有 效 自 动 屏蔽 了 宇宙 线 ， 在 流量 比较 低 的 
iim, (AMA LEAH aE EF IRAF 的 处 理 结果 。IRAF 需要 在 处 理 过 程 中 对 每 条 光谱 手动 屏蔽 宇宙 线 ， 相 
比 之 下 本 文 方法 节省 了 处 理 时 间 且 精度 更 高 。 
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图 8 结果 比较 。 蓝 色 曲 线 为 JRAF 处 理 结果 ， 黄 色 曲 线 为 本 软件 处 理 结果 

Fig.8 Comparison of results. The blue curve is the IRAF result, and the yellow curve is the result of this pipeline 
4 总 结 

本 文系 统 介绍 了 一 款 针对 2. 16 m 望远镜 BFOSC 光谱 仪 的 长 颖 光谱 数据 处 理 软件 ， 主 要 流程 包括 
生成 日 志 、 本 底 和 平 场 的 合并 与 改正 、 目 标 星 定位 、 视 场 畸 变 改 正 、 背 景 改 正 、 波 长 定 标 、 优 化 抽 
谱 。 本 文 的 数据 处 理 软件 使 用 方便 ， 相 比 TRAP 需要 手动 处 理 ， 本 软件 只 需要 用 户 输入 原始 数据 所 在 
的 路 径 即 可 ， 处 理 过 程 透明 可 视 ， 每 个 步骤 完成 都 会 生成 相应 的 图 片 并 保存 数据 以 便 重 复 使 用 ， 工 作 
效率 高 ， 能 批量 化 自动 处 理光 谱 数据 且 自 动 扣除 宇宙 线 数据 。 与 IRAP 处 理 结果 相 比 效果 更 好 ， 具 有 
较 高 的 实用 价值 。 

致谢 : 感谢 兴隆 2. 16 m 望远镜 全 体 工 作 人 员 的 支持 。 
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Abstract: In astronomical observation, the long-slit spectrometer is a kind of spectral observation equipment 


commonly used. The long-slit spectral mode of the Beijing Faint Object Spectrograph and Camera ( BFOSC ) 


mounted on the 2. 16 m telescope of the Xinglong Observatory, National Astronomical Observatories, Chinese 


Academy of Sciences, is widely used, but it is not equipped with a specific data processing software. After 


observation, users usually run IRAF software interactively to obtain the final data product. This paper 


introduces an automatic data processing software for BFOSC long slit spectrum. It includes bias and flat-field 


correction, background subtraction, spectrum extraction and wavelength calibration. We describe the 


processing method and the flow of the program, and present the results of processing with the observed data, 


and compare them with the IRAF results. Our software realizes automatic operation on the basis of completing 


all the necessary steps of spectrum processing, and adopts the optimal extraction method, which automatically 


deducts most cosmic rays and ensures the high precision of output results. 


Key words: long-slit spectrum; optimized extraction method; automatic pipeline 


